Polarization Evolution in Strong Magnetic Fields 
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ABSTRACT 

Extremely strong magnetic fields change the vacuum index of refraction. Although 
this polarization dependent effect is small for typical neutron stars, it is large enough to 
decouple the polarization states of photons traveling within the field. The photon states 
evolve adiabatically and follow the changing magnetic field direction. The combination 
of a rotating magnetosphere and a frequency dependent state decoupling predicts 
polarization phase lags between different wave bands, if the emission process takes place 
well within the light cylinder. This QED effect may allow observations to distinguish 
between different pulsar emission mechanisms and to reconstruct the structure of the 
magnetosphere. 



1. Introduction 

Understanding the structure of the magnetic fields surrounding neutron stars may provide a 
key in developing models for the radio and X-ray emission from pulsars, pulsar spindown, soft- 
gamma repeaters and the generation of the magnetic fields themselves. Although the magnetic 
field is instrumental in models of many phenomena associated with neutron stars, measuring its 
structure over a range of radii is problematic. Observations of the thermal emission from the surface 
may constrain the magnetic field geometry near the star, and the slowing of the pulsar's rotation 
may yield an estimate of the strength of the field near the speed-of-light cylinder. Connecting these 
regimes is difficult. 

The intense magnetic fields associated with neutron stars influence many physical processes - 



cooling ( Heyl fc Hernquist 1997c ; Shibanov et al. 1995 ), atmospheric emission ( Rajagopal, Romani 
fc Miller 1997 ; Pavlov et al. 1994 ) and the insulation of the core flHeyl fe Hernquist 1998b] |Heyl| 



fc Hernquist 1998a ; Schaaf 1990 ; Hernquist 1985 ). Even stronger fields such as thought to be 
found near AXPs and SGRs alter the propagation of light through the magnetosphere by way of 



quantum-electrodynamic (QED) processes and may further process the emergent radiation ( Heyl 
fc Hernquist 1997a| ; [Heyl k Hernquist 1997b| ; [Baring 1995| ; [Baring & Harding 1995| ; |Adler 197l[ ), 
and distort our view of the neutron star surface ( [Shaviv, Heyl Sz Lithwick 1999j ). 
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Even for weaker fields, QED renders the vacuum anisotropic. The speed at which light travels 
through the vacuum depends on its direction, polarization and the local strength of the magnetic 
field. Although for neutron stars with B < 10 14 G this effect is too weak to grossly affect images 
and light curves of these objects, it is strong enough to decouple the propagating modes through 
the pulsar magnetosphere. 

To lowest order in the ratio of the photon energy to the electron rest mass-energy, the 
index of refraction of a photon in a magnetic field is independent of frequency. Near the pair- 
production threshold, the photon propagation adiabatically merges with a postironium state 



( Shabad fc Usov 1986 ). However, well above the pair-production threshold for weak fields 



(m e c 2 <C E <C m e c 2 4.4x 10 13 G/B) the low energy results again provide a good approximation (|Tsai 



|fc Erber f975 ) In the context of the field surrounding a neutron star, only photons with sufficiently 



high wavenumbers (the optical and blueward) will travel through a portion of the rapidly weakening 
magnetic field without their two polarization modes mixing. 

At radio frequencies, the plasma surrounding the neutron star produces a similar effect 
(e.g. Cheng Sz Rudcrman 197£, Barnard 1986). One would expect that radio emission will be 



initially polarized according to the direction of the local magnetic field. When one observes a 
pulsar at a particular instant, one sees emission from regions with various magnetic field directions; 
therefore, one would expect the polarization to cancel out substantially. However, one finds that 
pulsars exhibit a significant linear polarization. As the polarized radiation travels from its source, 
its polarization direction changes as the local magnetic field direction changes. At a distance from 
the star that is large when compared to its radius, the local magnetic field in the plane of the sky is 
parallel across the observed portion of the star. Since the field changes gradually, the polarization 
modes are decoupled, and the disparate linear polarizations can add coherently. 

Previous authors have focussed on the propagation of polarized radio waves through the 
magnetosphere. At these frequencies, the vacuum polarization is safely neglected. Furthermore, 
they have assumed that the coupling of the two polarization modes occurs instantaneously. In 
this paper, we will treat the problem of vacuum polarization in particular and how the gradual 
coupling of the polarization modes affects the final polarization of the emergent radiation. At 
sufficiently high frequency the plasma only negligibly affects the radiation as it travels through 
the magnetosphere. The precise frequency at which the vacuum birefringence begins to dominate 
depends on the charge density of the magnetospheric plamsa. If we assume the Goldreich-Julian 



value ( Goldreich fc Julian 1969 ), we find that for 



the vacuum contribution to the birefringence dominates that of the plasma (the modes collapse 
where equality obtains ( [Meszaros I992j )). Our study of the vacuum-dominated regime, the optical 
and blueward, complements the work of |Cheng &: Ruderman (1979) and Barnard (1986) and will 



help to interpret observations of high-energy polarized radiation from neutron stars. 
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2. Small Amplitude Waves in the QED Vacuum 



] 



Kubo & Nagata (1981) (also [Kubo fc Nagata 1983] ) derive the equation of motion of polarization 
direciton on the Poincare sphere as light travels through an inhomogeneous birefringent medium. 
Since their results assume that the medium is polarized but not magnetized, we first extend their 
results to include magnetization. The general equations derived describe how polarized radiation 
travels through any birefringent medium in the limit of geometric optics. We then focus on the 
propagation of high frequency radiation through pulsar magnetospheres and how measurements of 
the polarization of this radiation can constrain both the structure of the magnetic field and the 
emission process of the radiation. This extended formalism is more than adequate to also describe 



the plasma induced birefringence of radio waves. This could be used to extend the works of Cheng 
& Ruder man (1979) and Barnard (1986) 



Kubo & Nagata (1983) find that the evolution of the polarization of a wave traveling through 
a birefringent and dichroic medium in the limit of geometric optics is given by 



ds 

dx 3 



— — = Axs+ Txs xs, 



(2) 



where X3 is the direction of propagation, s is the normalized Stokes vector ( [Jackson 1975| ), and 
Cl and T are the birefringent and dichroic vectors. The Stokes vector consists of the four Stokes 
parameters, So,Si,S2 and S3. The vector s consists of Si/So, S2/S0 and S3/S0. 

Waves traveling through a magnetized (or electrified) vacuum are best treated by linearizing the 
constitutive relations about the external field. |Heyl fc Hernquist (1997b) find that the permeability 
and dielectric tensors consist of an isotropic component plus an added component along the direction 
of the external field. Using the dielectric and permeability tensors of a magnetized vacuum in the 
formalism of Kubo fc Nagata (1983) yields 



kn 



2Jm£i + 1/2 (me f + fj,fEi) sin 2 6 



(fiiEf - fifEi) sin 2 9 



cos 2cj) 
- sin 2(f) 




(3) 



where 6 is the angle between the direction of propagation (k) and the external field and <p is the 
angle between the component of the external field perpendicular to k and the x— axis defined by 
the observer, Ui and [if are the isotropic and along the field components of the permeability tensor. 
£i and Ef are the equivalent dielectric tensors. 

For the case of vacuum QED with an external magnetic field to one loop, we find that the 
amplitude of ft is proportional to the difference between the indices of refraction for the two 
polarization states: 
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to lowest order in a, the fine-structure constant. £ = B/Bqeb (-Bqed ~ 4.4 x 10 G) and the 
functions Xq(x) and X\(x) are defined in Heyl fc Hernquist (1997a) . An external electric field 
yields similar results. However, in this case the vacuum is also dichroic so the vector T is nonzero. 



In the weak magnetic field limit (£ <C 0.5) we obtain, 

a 2 



n± 



Air 15 



¥ sin^ 6, 



and the strong field limit (£ 3> 0.5) yields 



n± 



nil 



a 2 
4^3 



£sin 2 0. 



(6) 



(7) 



2.1. Exact Solutions 

Kubo & Nagata (1983) found that Equation Q yields exact solutions for restricted values of Cl 
and T. In the case of a uniformly magnetized vacuum, T = and Ct is constant. In this case, the 
polarization vector s traces a circle on the Poincare sphere about the vector 12 at a rate of \€l\. 

I . I 

Kubo & Nagata (1981) examine the case where ft is constant in magnitude but satisfies, 

|^ = r x Ci (8) 

dx 3 

where T is a constant. In this case Cl rotates about T at a rate of |T|. If we follow the equations 
in a rotating coordinate system such that in it ft' is constant, we find that s' satifies the following 
equation 

|^ = S\ ff x s' (9) 

where 

{i' cS = ti'-r. (io) 

This equation holds even if T is not constant. However, if T is constant, we obtain a new exact 
solution where s circles a guiding center displaced from ft which in turn rotates about T. 

If we take s || f2 initially and T _!_ Cl, we find that s develops a component perpendicular to f2. 
In the case of vacuum QED, ft lies in the 1 — 2— plane. If the magnetic field rotates uniformly in 
the plane transverse to the wave, we find that s will leave the 1 — 2— plane and a circularly polarized 
component will develop. 



2.2. Adiabatic Approximation 



If the parameters describing the motion of a system vary slowly compared to the characteristic 
frequency of the system, the system evolves adiabatically such that at a given time it executes a 
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motion given by the instantaneous values of the parameters as if they were static. The exact solution 
given by Equation |9| has this feature if \ ft\ S> |T|. If this limit applies, Cl e g is nearly parallel to Cl' 
and s circles the instantaneous guiding center Cl. Furthermore, if s is initially parallel or antiparallel 
to Cl, it will remain so (corresponding to polarizations that are parallel or perpendicular to the 
direction of the magnetic field). That is, the polarization modes are decoupled, and the polarization 
direction follows the direction of the magnetic field. 

We can also build an adiabatic approximation onto the exact solution of Equation ^ by allowing 
the magnitude of Cl to vary. We take a wave polarized parallel to the initial value of Cl. As Cl 
rotates about T and decreases in magnitude, s rotates about f2 e g- and follows the instantaneous 
direction of f2 e ff as long as 

' 1 d\Cl\ 



ft 



\Cl\ dx 3 



» 1. 



(11) 



Numerical integration of Equation ^ for \Cl\ = Ax^ G and T = — 1/10x3 (Cl rotates by one radian 
after the photon has travelled 10 units of distance) shows that the polarization follows the analytic 
solution described in the previous paragraph for 



1 d\Cl\ 
\Cl\ dx 3 



£0.5 



(12) 



and then freezes for values of A ranging from 10 to 10 8 . Figure |] depicts both the numerical results 
and the analytic approximation. 

Generically, if the polarization modes of the medium are linear, and the wave begins in one of 
the modes, the final direction of the polarization projected into the S± — S2— plane will depend on 
how much the field direction has changed by the time of mode coupling. The circular component of 
the polarization depends on how fast the field direction is changing at the time of mode coupling. 

By applying these results to a magnetic dipole in vacua, we find that an outgoing photon's 
polarization will follow the analytic solution until 

SdlnAnS 1 / 3 
2 dlnS An 



k Q R B 



1/3 



(13) 



after which its polarization in the observer's system freezes. The left side of the equality describes 
the environment at the point when the polarization freezes; the right side depends on the point of 
emission of the photon. Figure ^ shows the magnetic field strength at coupling as a function of 
photon energy for typical pulsar parameters. 

If the modes begin to mix where B <C -Bqed (this is appropriate for all but the lowest energy 
photons near the most strongly magnetized neutron stars), we obtain 
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cm 



If this ratio is less than unity for a photon, the photon will travel with its polarization modes 
decoupled during some portion of its journey away from the star. That is, if 

£>M9xl0 -3 eV (^)- 2 (_|L)-\ (16) 

the vacuum will decouple the polarization modes. For parameters typical to neutron stars, the 
polarization modes for all photons with A <J 800 ^m will be decoupled at least as they begin their 
journey from the vicinity of the neutron star. Although for photons of such low energy, the mode 
coupling induced by the plasma is likely to dominate (e.g. Cheng fc Ruderman 1979] ) . 



These effects may also be important for photons travelling through the magnetospheres of 
strongly magnetized white dwarfs. The most strongly magnetized white dwarfs have B ~ 10 9 G 
and Rq ~ 10 9 cm; in the magnetospheres of these stars, the polarization modes of photons more 
energetic than 1.5 eV will be decoupled by the vacuum birefringence. 



3. Rotating Pulsar Magnetospheres 

The photons that we observe from neutron stars have to travel through a portion of the neutron 
star's magnetosphere before reaching us. If they travel through a region where Equation [ll] holds, 
the photon's polarization relative the observer's axes will change as the magnetic field direction 
changes. If we examine a photon which leaves the surface of the star, once it has travelled a 
distance of several radii it will be travelling approximately radially. If the magnetic moment of 
the star remains fixed during the journey and the modes remained decoupled, the polarization 
directions of all the photons travelling in a particular direction and polarization mode will align 
with each other. 

The plasma in the vicinity of a neutron star is polarized by passing waves of sufficiently low 
frequencies. |Cheng fc Ruderman (1979) used this effect to explain the high linear polarization 



of radio emission from pulsars even though one would expect the polarization modes where the 
radiation is produced to vary from place to place. The same effect will operate at high frequencies 
through the vacuum polarization of QED. 

High frequency photons travelling through a pulsar magnetosphere may travel a large distance 
before the modes couple. During this travel, the magnetic field of the pulsar may rotate a significant 
amount, thereby changing the polarization modes in the observer's frame as the photon travels. 
Since high energy photons travel further than less energetic ones, the higher energy photons will 
have their polarizations dragged further by the rotating magnetic field. 
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3.1. Analytic Treatment 

If we assume that the rotation rate of the projection of the magnetic field onto the plane 



transverse to a photon's propagation changes sufficiently slowly, we can apply the results of § ^2 
to calculate the final polarization of the photon. What is required is the direction of magnetic field 
at the point of recoupling and its rate of change. 

This analytic treatment will be restricted to photons travelling in the radial direction. For 
a significant phase lag to develop, the polarization modes must remain decoupled until well away 
from the star. Here outgoing photons emitted from near the stellar surface will have approximately 
radial trajectories. Furthermore, at these distances both gravitational light bending and magnetic 
lensing ( [Shaviv, Heyl fc Lithwick 1999| ) may be neglected. 

The photon's trajectory will be characterized by the angles £ and 9 the longitude and colatitude 
where the photon left the star. We shall take w, the phase of the star's rotation, to be zero when 
the photon leaves the surface. The magnetic pole lies at zero longitude and colatitude \P. As the 
photon travels away from the neutron star, the phase w increases and the local magnetic field 
direction changes. The criterion for recoupling in the weak-field limit when travelling through a 
dipole field is 




0, equator \ D 6 , s ^ a 

— K k — — 

-DQED / r v 



dsin a 1 6 



dr sin 2 a 



\ < 17 > 



In this equation, the angle a designates the magnetic colatitude of the photon which changes as 
the star rotates underneath it. It satisfies the following expression 

cos a = sin 9 cos(£ + w) sin ^ + cos ^> cos 9. (18) 

The second important angle is the angle between the observer's axes and the local magnetic field 
direction in the plane of the sky. For a magnetic dipole, the projection of the local magnetic 
direction and the projection of the magnetic moment (m) onto the tangential plane are colinear. 
Unless the line of sight is directed down the rotation axis, we can use the projection of the rotation 
axis (z) as one of the reference axes for measuring the polarization of the radiation. These projected 
vectors are defined by 

z p = z — f cos 9 (19) 
nip = m - f cos a. (20) 

The travelling photon does not distinguish between a magnetic field pointing in one direction or in 
the exactly opposite direction (see Equation^); therefore, we only need to know the angle between 
z p and m p modulo n. Calculating the cross product suffices to give the magnitude of the angle, 



which yields 



(21) 



sin(f + w) . . 

tan0= -, ^4 — n- ( 22 ) 

cos 9 cos(£ + w) — cot W sin 9 
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T in this case is given by 



deb 2n ~ 
dm cF 



(23) 



So the final polarization of the photon in the observers frame is given by Equation 1C evaluated at 
the moment of recoupling. 

The calculation of the polarization evolution in the analytic treatment proceeds as follows, 



1. Choose the colatitude of the observer (0), the longitude where the photon is emitted (£), a 
reference frequency and magnetic field strength. 



2. Given the period of the pulsar (P), a photon's radius is given by 

r = —w + R . 
Zn 



(24) 



For a given value of w, the phase, solve for Po,equator^o assuming that Equation |17] is satisfied, 
and calculate the magnitude of fl at recoupling. 



3. Substitute this value of w into Equation 22 to calculate the direction of ft at recoupling. 



4. Finally, calculate T at recoupling using Equation 23. The final polarization according to this 
analytic adiabatic approximation lies along 



cff 



n-r. 



(25) 



5. This procedure can be repeated for other values of the longitude (£) to generate a light curve. 



In the case where the line of sight is aligned with the rotation axis, many of these geometric 



considerations vanish, and the problem reduces identically to that treated in § 2.2. 



In this situation we can derive expressions for both the final position angle, 
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and the circular component, 
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(26) 
(27) 
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If | S3 1 / 1 s I <C 1, the following approximation holds, 



For a more general pulsar geometry we can still use Equation 26 to characterize how radiation 
at various frequencies is polarized. Specifically, the position angle can be converted to a time lead 
for the polarization angle of high-frequency relative to low-frequency radiation. This time lead 
given by P0/(2tt) is independent of the period of the pulsar. Figure || plots the known radio 
pulsars ( Taylor, Manchester fe Lyne 1993; ) and the approximate phase lead (r coup i e /ri c ) expected 



at 5.2 keV relative to zero-energy photons. If r coup i e is a large fraction of the radius of the light 
cylinder, it may be possible to study the geometry of the magnetic field at large distances from the 
pulsar. For the Crab pulsar the radius at coupling for 5.2 keV is about 16% of the light-cylinder 
radius. The millisecond pulsars, PSR J1939+2134 and PSR J1824-2452, have r couple of over 9% of 
n c at 5.2 keV. 



3.2. Numerical Treatment 

The analytic techniques outlined in the previous subsection provide insights to interpret and 
extend the results from a direct numerical integration of Equation ||. Nevertheless, a numerical 
integration is unavoidable if we wish to calculate the phase shifts under general conditions - when 
the observer's inclination angle is nontrivial and the field is not strictly dipolar. 

The numerical integration assumes two basic assumptions: 

1. We assume that either 

(a) The magnetosphere is co-rotating with the NS. This implies that the local magnetic field 
is at any given moment aligned with the NS, irrespective of the distance. (A deviation 
from this behavior will indicate where the co-rotating magnetosphere ends), or 

(b) The magnetic field is given by that of a magnetized conducting sphere rotating in vacua 



(Deutsch 1955, Barnard 1986) 



2. The photon is traveling radially. That is to say, the emission process takes place well within 
the region where the polarization states couple. 

With these assumptions considered, the integration is achieved in a straightforward manner. A 
photon is followed from an initial radius i?o to a radius rfj na i that is beyond the recoupling distance 
using the equation: 

ds 

- = flxs. (31) 
or 

Each radial step is integrated forward using the fourth order Runge-Kutta algorithm. To do so, one 



has to calculate the birefringent vector fl using Equation 22 for the angle of the magnetic field and 
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Equation || for the strength of the birefringence. To calculate this vector for the Deutsch (1955) 



fields, we use Eq. 7 of Barnard (1986). The rotation phase of the star is related to the radial 
coordinate through Equation 24. 



3.2.1. The Crab pulsar 

The Crab pulsar is among the most thoroughly studied astronomical objects. Its fast period 
and moderately strong magnetic field make it an ideal example for this process. We take 6 = 54° 
and ^ = 64° which reproduces the magnitude of the polarization swing for the pulse and interpulse 
(Smith et aL 1988), and the field strength at the magnetic equator to be 1.9 x 10 12 G (Taylor 



Manchester &: Lyne 1993] ) . The results of this calculation are depicted in Figure || and Figure ^. If 



the field is indeed a corotating dipole, the delayed coupling of the polarization modes results in the 
polarization of high frequency radiation leading that of lower energies. The lead time estimate from 



Equation |26j works well away from the pulse or interpulse where it slightly over or underestimates 
the lead time. However, if the neutron star is surrounded by a Deutsch field, this lead is retarded 
slightly. 

Since the coupling of the modes occurs gradually, a significant circular component will develop 
when the modes are weakly coupled. If the radiation initially has a circular component, this initial 
component is washed out when one measures the circular polarization over a finite bandpass, leaving 
only the value of S3 produced by the coupling. Both the Deutsch fields and the corotating dipole 
fields result in a significant amount of circular polarization (up to 30% of the total polarization 
for the corotating dipole at 5.2 keV or 7% of it at 4.5 eV). The circular polarization generated by 
the rotating Deutsch field peaks at about 65% of the total polarization at the highest frequency 
studied. 



3.2.2. RX J0720. 4-3125 
RX J0720. 4-3125 is an isolated neutron star candidate suspected to be close to the Earth. 



Heyl fc Hernquist (1998c) from cooling arguments and its current period of 8.391 s estimate that 
it has a magnetic field of B ~ 10 14 G. From Equation 26 we estimate the time lead for a given 



frequency to be approximately five times larger than for the Crab pulsar. As we see from Figure)^, 
since the period of RX J0720. 4-3125 is much larger than that of the Crab, the expected change in 
position angle over the larger lead time is very small. Furthermore, since the coupling of the modes 
occurs well within the light cylinder, the difference between the Deutsch and the corotating dipole 
model is negligible. Since the circular polarization that develops during coupling is proportional 
to the rotation frequency of the pulsar, very little circular polarization is evident in the emergent 
radiation - S3 peaks at 2.7% of the total polarization, coincident with the pulse and interpulse. 
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4. Discussion 

The vacuum polarization induced by quantum electrodynamics decouples the polarization 
modes in the strong magnetic field surrounding neutron stars. As radiation travels away from the 
neutron star, the modes of low- frequency radiation couple first; consequently, the position angle of 
low frequency emission will lag behind that of higher frequency radiation emitted coincidently and 
in the same mode. Puslars often exhibit mode switching as a function of phase and frequency. An 
inelegant way to avoid confusion is to examine the position angle modulo ir/2. By measuring the 
position angle as a function of frequency, one can determine not only the location of the emission 
but also probe the structure of the magnetic field near the light cylinder using sufficiently high 
frequency radiation. 

Previous authors have calculated a similar effect where the plasma decouples the modes at radio 
frequencies, assuming that the coupling occurs instantly ( Cheng &: Ruderman 1979| ; Barnard 1986). 



To relax this assumption, we use the equations describing the propagation of polarized radiation 
through a polarized and magnetized medium, under the assumptions of geometric optics. These 



results extend those of [Kubo fc Nagata (1981)| to the case of a magnetized and polarized medium. 
Linearly polarized radiation will develop a circular component during the gradual coupling of the 
modes. 

In the case of QED in a vacuum, the propagating modes are purely linear; therefore, the 
decoupling will wash out any initial circular component averaged over a finite bandpass. The only 
circular component present in the outgoing radiation is produced during coupling and measures the 
ratio of the coupling radius to the radius of the light cylinder. This circular component provides 
an independent validation of the effect. 

Since high frequency radiation is generally measured incoherently, determining the Stokes 
parameters becomes more difficult with increasing photon energy. Blueward of the ultraviolet, 
measuring circular polarization is impractical. The Spectrum- X-Gamma mission will carry a stellar 
x-ray polarimeter (SXRP) which will be especially sensitive at 2.6 keV and 5.2 keV, the first and 
second-order Bragg reflections off of graphite. Future instruments may be able to measure both 
linear and circular polarization at yet higher energies. 

From a theoretical point of view, the study of the propagating modes of a strongly magnetized 
vacuum is simpler than analysing those of a strongly magnetized plasma, especially when the 
properties of the plasma are not well known. However, given the magnetic permeability and the 
dielectric permittivity of the circumpulsar plasma, these results can be applied to study plasma- 
induced decoupling on the properties of polarized radio emission from pulsars. We leave the details 
of this process for a subsequent paper; however, our general arguments apply here as well. 

Both the plasma-induced and vacuum-induced decoupling can be used to probe the structure of 
pulsar magnetospheres. A measurement of the circular polarization produced by the mode coupling 
yields an estimate of the radius of the coupling. The frequency corresponding to a given degree 
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of circular polarization measures the radial gradient of third power of the difference between the 



indices of refraction of the two modes. This results from the approximate inequality Equation 12 



Like the arguments used in § 2.2, it is generic to waves travelling through a birefringent medium 
in the limit of geometric optics. 

Unlike the phase lead effect found by previous authors for the case of the plasma alone, the 
circular polarization provides a local probe of the pulsar magnetosphere without having to compare 
observations over a wide range of energy. As long as the radiation is produced in a region where 
the modes are decoupled, the interpretation of the circular polarization is straightforward. 

Measurements of the position angle of the radiation as a function of phase and frequency can 
also yield the same information, since the induced circular polarization is simply related to the 
phase lead. Furthermore, measuring all four Stokes parameters can elucidate the location of the 
emission process since below a critical frequency the modes will be coupled at the point of emission 
and the observed position angle will be constant with frequency. The circular polarization of the 
first modes which do travel through a decoupled region yields the radius of the emission. The 
details of the interpretation will only depend weakly on the medium which decouples the modes - 
plasma or vacuum. 

We have derived equations which describe the propagation of polarized radiation through a 
magnetized and polarized medium. By applying these results to the strongly magnetized regions 
surrounding neutron stars, we have found several observational probes of pulsar magnetospheres 
and emission. 
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Fig. 1. — Numerical solution and analytic approximation. The figures depict the exact numerical 
solution (solid lines) and the analytic solution based on the adiabatic approximation (dashed lines). 
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Fig. 2. — The dipole magnetic field strength at coupling. The strength of the magnetic field at 
coupling depends on the magnetic moment of the star and the energy of the photon. The modes 
of higher energy photons couple later (further from the star) . 
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Fig. 3. — The observed polarization as a function of frequency for a line of sight parallel to the 
rotation axis in the presense of a rotating dipole field. The panel depicts the polarization angle 
(solid line) and circular component (dashed line, normalized to |s|) as a function of wavelength. 
The lines from top to bottom trace the solutions for P = 1, 10, and 100 ms respectively. 
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Fig. 4. — The ratio of the coupling radius at 5.2 keV to the radius of the light cylinder for the 
known radio pulsars. Only pulsars with r coup ii ng /ri c > 0.01 are plotted. 
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Fig. 5. — The evolution of the Stokes parameters on the Poincare sphere for photons leaving the 
surface of the Crab pulsar at the moment the observer crosses the Prime Meridian. The bold 
lines depict the corotating dipole, and the light ones the Deutsch field. The energies plotted are 
5.2 keV and 5.2 MeV. The bold line in the 1-2 plane is the birefringent vector f2. The photons' 
polarization vectors and the birefringent vector all start from the -1 direction. The lower energy 
photons follow f2 for a shorter time than do the higher energy ones. Since the mode recoupling is 
not instantaneous, a circular component (S3) is generated. 
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Fig. 6. — The position angle as a function of photon energy and time for our model of the Crab 
pulsar. From left to right, the energies are 5.2 Mev (dot-dashed), 5.2 keV (long-dashed), 2.6 keV 
(short-dashed), 4.48 eV (dotted) and zero energy (solid). The bold lines follow the corotating dipole 
model, and the light lines follow the Deutsch model. 
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Fig. 7. — The difference between the final position angle of a photon emitted at a particular time 
and energy and the position angle of a zero energy photon emitted at the same time. The model 
is for the Crab pulsar surrounded with a Deutsch field. 
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Fig. 8. — The circular component of the radiation as a function of photon energy and time for 
the Crab pulsar, assuming the initial polarization is complete (i.e., |s| = 1; otherwise, it should be 
normalized to its smaller value). The energies are 5.2 MeV (dot-dashed), 5.2 keV (long-dashed), 
2.6 keV (short-dashed), 4.48 eV (dotted) and zero energy (solid). The bold lines follow the 
corotating dipole model, and the light lines follow the Deutsch model. 
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Fig. 9. — The position angle as a function of photon energy and time for our model of RXJ 0720.4- 
3125. We have taken -B ,equator = 10 14 G, = 85° and # = 90°. From left to right, the energies are 
5.2 Mev, 5.2 keV, 2.6 keV, 4.48 eV and zero energy. 
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Fig. 10. — A schematic depiction of the expected phase lag at a given rotation phase of the NS. 
Lower frequency photons follow the changing magnetic field only over a short distance from the 
neutron star, and thus the modes couple when the local direction of the magnetic field has not 
changed much since emission. The polarization modes of higher energy photons couple further 
away and thus follow its changing direction further, creating a phase lead between the photon and 
the NS dipole direction. The magnetosphere rotates as a solid body only up to the co-rotation 
radius (which must be less than the radius of the light cylinder); consequently, a photon's phase 
lead will not be increased anymore once it begins traveling through the retarded magnetic field. 
If the emission process takes place at a finite radius, the lower energy photons could already be 
emitted outside the decoupled region. In such a case, they would suffer no phase lag. 



